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Finally,  the  extension  of  the  MFE  formulation  to  three-coordinate 
parabolic  systems  and  for  the  transonic  small  disturbance  equations 
is  briefly  described. 
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1 . Introduction 

An  integral  procedure  using  spline  polynomials  is  described 
for  the  two  dimensional  boundary  layer  equations.  This  is  a modi- 
fied finite-element  (MFE)  formulation,  wherein  each  term  in  the 
equations,  rather  than  each  independent  variable,  is  approximated 
with  a spline  curve  fit.  Therefore,  this  is  not  a true  finite- 
element  or  Galerkin  method  and  the  conventional  spline  relationships 
between  functional  and  derivation  values  still  apply.  The  only 
difference  between  the  present  integral  formulation  and  our  earlier 
differential  collocation  procedures  is  in  the  treatment  of 

the  governing  differential  equations.  The  differential  methods  are 
more  suited  to  non-conservation  equations;  the  present  integral 
formulation  is  more  desirable  for  conservation  or  divergence  form  of 
the  equations . 

For  the  boundary  layer  equations  in  divergence  form  it  is  shown 

here  that  the  spline  MFE  approach  is  equivalent  to  the  class  of  two- 

(4) 

point  methods  proposed  by  Keller  . His  second-order  development 
commonly  termed  the  Keller  Box  Scheme*  (KBS)  has  been  widely  used 
for  boundary  layer  investigations.  It  has  been  suggested  that  due  to 
the  two-point  centered-dif ference  character  of  this  formulation,  it 
is  second-order  accurate,  even  for  non-uniform  grids.  In  the  present 
paper  it  is  shown  that  when  the  KES  difference  system  is  reduced  to 
an  equivalent  three-point  system,  the  governing  equations  are  in  fact 


*The  two-point  KBS  development  reproduced  herein  is  equivalent  to  the 
original  half-point  difference  version  of  the  KBS  only  for  equations, 
in  divergence  form.  This  point  is  not  clearly  expressed  by  Wornom 
in  his  applications  of  the  KBS,  see  Section  7. 
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only  satisfied  to  first-order  in  the  usual  finite-difference  sense; 

1. e.,  when  the  truncation  error  is  evaluated  at  the  mesh  point,  . 

It  is  further  shown  that,  for  a model  equation,  a simpler  procedure 
for  obtaining  the  tridiagonal  form  of  the  KBS  equations  can  be  de- 
rived from  the  MFE  integral  method,  or  alternatively  a three-point 
weighted  averaging  procedure.  These  lead  directly  to  a single  tri- 
diagonal  structure,  for  a second-order  differential  equation,  in 
place  of  the  2x2  block-tridiagonal  system  resulting  from  the  usual 
KBS  formulation.  The  MFE  approach  is  extended  to  higher  order  by 
using  cubic  splines  and  it  is  shown  that  this  leads  to  the  fourth- 
order  version  of  the:  KBS  procedure  recently  applied  by  Wornom^^^  . 

Boundary  layer  solutions  using  conventional  second-order  finite- 

difference  collocation,  the  second  and  fourth  order  KBS,  and  fourth- 

(3,9^ 

order  spline  collocation  ' or  MFE  methods  are  compared-  Conservation 
and  non-conservation  form  are  considered.  Finally,  the  extension 
of  the  MFE  formulation  to  three-coordinate  parabolic  systems  and  to  the 
transonic  small  disturbance  equations  is  briefly  described.  Solutions 
with  this  development  have  not  yet  been  obtained. 

2 . Keller  Box  Scheme;  Conservation  Equations 

Let  us  apply  the  KBS  method  to  the  following  ordinary  differen- 
tial equation  written  in  divergence  form  and  with  a source  term  F: 

(au,)  ^ = F + (bu)  ^ ;a  = a(u^,u;'T');F  = F(u;“)  (1) 

b = b (u;  P)  . 


The  KBS  met.hod  considers  (1)  as  two  first-order  equations, 
u ^ = m 

(am)  ^ = F + (bu) 


so  that, 


(2a) 


2 


r 


:2b) 


The  KBS  finite  difference  anajog  of  (2a)  is: 

h . 

u . -u . T = ^ (m.+m.  t) 

: D-1  2 j j-1 

h . 

(am)  (am)  . , = (F  .+F  . , ) + (bu)  (bu)  . ^ 

J 3-i  2 3 3-1  'j  j-1 

where  h . = . 

1 1 D-1 

Note  that  the  right-hand  side  of  (2b)  ic  simply  the  trapezoidal  rule 

j 2 

applied  to  f Fdr  . The  error  in  (2b)  at  (j-^)  is  formally  0(h  ). 

■ ”1-1 

This  discretization  can  be  interpreted  as  the  application  of  a linear 
polynomial  for  the  integrand.  Other  polynomial  forms  will  lead  to 
different  formulations  and  can  provide  higher-order  extensions  of  the 
KBS  method.  In  a subsequent  section  we  will  derive  one  such  scheme 
and  show  its  equivalence  to  the  fourth-order  KBS  applied  by  Wornom. 

2 . 1 Reduction  of  KBS  to  Tridiagonal  Form. 

Let  us  eliminate  the  intermediate  function  m^  from  (2b).  Solving 
for  tftj  and  m^  we  find  on 


(3a) 


+ (bu)  j- (bu) 


(a .+a  . , )m . , 
3 3-1'  3-1 


(bu)  (bUj_j^; 


(3b) 


Similar  expressions  for  m.  and  m.  , are  found  on 
^ 3 3+1 

crementing  j . After  some  elimination  we  find 


by  in- 


2a • . , 2a  . T 

J±i  _ —2rl 


•^j+1 


-(bu)  ,^j+(bu)  = (a.^^-a^_j)m. 


2 

(4a) 


3 


with 


= r . b = constant  and  a^  = constant  = 1,  we  see  that 


"(l+7)h. 


2 2at^) 


It  should  be  noted  that  for  non-uniform  grids  {7^1),  both  u^^  and 

u^  in  (1)  are  represented  by  first-order  accurate  finite-difference 

discretizations.  If  the  truncation  error  is  formally  evaluated,  at  T., 

3 

this  three-point  analog  of  the  KBS  method  is  in  fact  only  "first-order 
accurate".  The  most  interesting  aspect  of  equation  (4b)  is  the 
weighted  average  of  the  source  term  F;  this  implies  a certain  amount 
of  smoothing  not  found  in  a differential  collocation  method.  This 
three-point  reduction  has  previously  been  used  by  Ackerberg  ^ ' and 
others.  It  has  not  been  extensively  applied  and  appears  to  be  more 
efficient  than  the  2x2  KBS  solution  procedure.  For  variable  a ^ , as 
in  the  case  of  turbulent  flows,  m^  appears  explicitly  in  (4a)  and 
can  be  eliminated  with  (3a)  . The  tridiagonal  system  is  more  complex. 

3 . Modified  Finite-Element  Formulation  (MFE) 

Equation  (4b)  can  be  derived  directly  by  integrating  equation  (1) 
over  two  adjacent  intervals  so  that, 

, j '3+1 

I [u^^- (bu)  ^-F]dTi  + f Cu^^^- (bu)  ^-F]dT=0  ; m = u^  (5) 


The  variables  u ( *;)  and  F (?)  are  then  approximated  by  a polynomial  over 
each  interval  C j _ » “n j 1 » j=l»  .../N.  For  example,  with  a quadratic 
polynomial  approximation  for  u over  each  of  the  intervals,  and  a 
linear  polynomial  consistent  with  the  trapezoidal  rule  of  integration, 
for  the  source  term  F*,  we  obtain  on  [ ^ ?^]: 


*A  lower  order  polynomial  for  F is  consistent  with  a quadratic  approxi- 
mation for  derivatives.  The  quadratic  polynomial  in  (6)  has  been 
described  earlier  in  reference  (2)  and  is  designated  as  S(?;2,l). 


u(Tl)  = Ujt+Uj_j^  (1-t)  + (Uj-Uj_^-hjmj)  t(l-t) 
F(T1)  = Fjt+Fj_j_(l-t) 


where  t = 


All  fiinctional  values  appearing  in  (5)  are  assumed  to  be  individually 
continuous,  i.e.,  u,F,m. 

Completing  the  integration  in  equation  (5) , we  obtain 

h . 

^"'j+l-"'j-l^-^(''j+l-’"j-l^+  2^  C7F.^^+(l+a)Fj+F._^:=0  (7a) 

Differentiating  (6)  we  recover  the  expression  in  (2b) 


u . -u  . , 
3 3-1 


2 


The  expression  (7b)  is  also  obtained  by  imposing  continuity  of  m^ 
Eliminating  m^  from  (7a, 7b)  we  find 


u . -u  . , 
1+1  1-1 


(7  (1+a)  h . 


2 ^ (l+a)h 


- {cFj_^j_+(l+a)Fj+Fj_3_)/(2  (l+a); 

This  equation  is  exactly  the  same  as  (4b)  and  as  mentioned  previously 

has  the  advantage  of  being  of  a single  tridiagonal  form.  This  system 

can  more  easily  be  solved  by  the  usual  two  pass  algorithm  than  can 

the  2x2  block  inversion  for  the  KBS  method  and  provides  identical 

solutions.  For  F-0,  (7c)  is  equivalent  to  the  difference  form  of  the 

(7) 

momentum  equation  in  the  Davis  Coupled  Scheme 

This  modified  finite-element  formulation  differs  from  the  usual 
finite-element  method  in  that  the  weighting  functions  are  unity  and 
olvnomials  are  prescribed  for  each  term  in  the  equation  rather  than 


for  each  variable.  This  leads  to  a considerable  simplification  and 
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provides  consistent  error  estimates  for  the  method.  Finally,  this 
procedure  also  applies  for  a variable  coefficients  aj,bj;  in  which 
case  (4a)  is  recovered;  it  can  be  made  more  accurate  by  consider- 
ing higher-order  polynomials.  This  is  shown  in  Section  6.  The 
polynomial  for  u ( ?)  in  (6)  leads  to  the  usual  finite-difference  re- 
lationships if  instead  of  (7a),  we  require  the  continuity  of 
This  leads  to  3(772,0)  of  reference  (2).  The  MFE  approach  does  not 
require  this  additional  continuity  condition. 

4 . Weighted  Averaging 

The  usual  three-point  differencing  for  first  and  second  deri- 
vatives is  given  by: 


(8a) 


and 


7 ( l+g 


(8b) 


In  terms  of  the  mesh  half-points,  these  difference  formulas  can  be 
interpreted  by  the  following  weighted  averages: 


u_  = m . 


m . , I +^m  . 1 
(1+r) 


(9a) 


where 


u,  = 

3 l+(7 


u . -u . 


U . -U  . 1 


(9b) 


(10) 


6 


assumed  in  lieu  of  (9),  we  have  the  following  formulas  for  : 


m . = u _ = 


M.  = 

: irn 


(1+ct) 


l+a 


(12) 


With  the  same  averaging  procedure  for  all  source  terms,  and  the 
additional  conditions  (10)  and  (11)  , we  recover  the  tridiagonal  form 
of  the  KBS  or  modified  finite-element  method,  equation  (4) . It  is 
interesting  that  as  with  the  modified  finite  element  approach,  the 

! 

I weighted  averaging  here  is  applied  to  the  governing  equation  and  not 

the  individual  variable. 

If  the  expression  (9)  for  is  left  in  terms  of  m^  ( from  (11)), 
the  2x2  (u,m)j  system  is  truly  second-order  accurate  for  both  m^  and 
M j ; the  trapezoidal  averaging  is  only  first-order  accurate.  Solutions 
using  this  formulation  have  never  been  considered. 

As  a final  note,  we  have  conducted  some  sample  boundary  layer 
calculations,  both  laminar  and  turbulent,  with  and  without  .a  pressure 
gradient.  These  results  are  shown  in  Section  7.  It  is  significant 

1 

) that  neither  the  downstream  weighted  formulation  (KBS  or  MFE) , nor  the 

I upstream  weighted  formulation  (usual  three-point  differencing)  is 

superior  throughout.  The  preferred  method  depends  on  both  the  flow 
profile,  and  significantly,  the  finite-difference  grid.  Therefore, 
once  again  no  single  formulation  can  be  interpreted  as  a panacea,  even 

j for  a limited  class  of  flows,  i.e.,  simple  boundary  layers. 
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5 . Fourth-Order  Keller  Box  Scheme 

Recently  Wornotn  has  used  the  fourth-order  extension  of  the 

(4) 

KBS  proposed  by  Keller  . The  extension  of  the  usual  trapezoidal 
rule  to  fourth-order  gives  the  following  two-point  Taylor  series 


expansion : 


When  applied  to  equations  (2a) , we  find 


U . -U  . , = (m  .+!D  . , ) 

3 3-1  2 ' :-l' 


12 


(am) j- (am) j_^=(bu) j- (bu) (Fj+Fj_^; 


" i ^ ^ I 

rt  J 

The  primes  denote  differentiation  with  respect  to  T.  The  expression 
(14a)  is  equivalent  to  one  of  the  spline  collocation  formulas,  see 
references  (1-3)  and  Section  6.  As  applied  by  Wornom  , the 
system  (14)  is  closed  by  evaluating  Mj  directly  from  the  governing 
equations.  This  results  in  a 2x2  block  for  (u.m)^.  An  alternate 
procedure  shown  in  the  next  section  will  be  to  apply  the  spline  rela- 
tionships between  m^  and  . 

If  we  increment  j to  j+1  and  add  the  resulting  equations  to  (14) , 
we  obtain  a three  point  version  of  this  scheme.  The  resulting  equa- 


tions are  as  follows: 


h . 

-J.  r. 


8 


(16) 


h. 


[anVj+i+U.+  c)m.+m._^] 


-h.' 


1 


12  '3+1  -:■  -3  -3- 

5.1  Keller  Box  Schemes;  Non-Divergence  Equations 

If  we  neglect  the  higher  order  terms  in  the  two-point 
formula  (13),  we  have 


h.  , , 


'3  -3-1  2 --3  -3 

As  noted  earlier,  this  is  equivalent  to  the  application  of  the 
trapezoidal  rule  and  when  applied  to  equation  (2a)  leads  to  the 
second-order  KBS.  However,  the  application  of  the  above  algorithm 


(17) 


to  the  non-divergence  form  of  the  equations  does  not  lead  to  the 

(8) 


original  KBS  method  . Equations  (2a) , with  a = 1,  can  be  re- 
written in  non-divergence  form  as 


= m 
T| 


m_  = bm+F  where  F = F+b_u. 
T1  T| 


(18a) 

(18b) 


Application  of  (17)  to  this  set  of  equations  leads  to 


h. 


and 


u.-u.  , = (m.+m. 

33-1  2^3  3-1^ 

h. 


(19a) 


mj-mj-i 


(b  . m . +b  . , m . , +F  . +F  . , ) 
2 33  3-1  3-1  3 3-1 


(19b) 


Using  half-point  differencing,  as  is  done  in  the  original  KBS  method, 
the  equation  (18b)  becomes 


h. 


b .+b  . 


m 


(19c) 


It  is  easily  seen  that  equations  (19b)  and  (19c)  are  identical  only 
for  constant  b,  i.e.,  divergence  form.  The  calculations  of  section  7 


will  show  that  in  certain  cases  there  can  be  a significant  increase 
in  accuracy  when  applying  (19c)  in  place  of  (19b)  . The  reduction 
to  tridiagonal  form  follows  the  same  course  and  for  variable  grids 
the  system  remains  first-order  accurate.  For  divergence  form 
equations  the  half-point  differences  and  two-point  series  formulas 
are  identical.  The  accuracy  of  the  fourth-order  two-point  formula 
(13)  for  non-divergence  equations  has  not  been  established  as 
Wornom  only  considers  the  divergence  form  equations  characterized 

by  (14b).  If  non-divergence  form  is  considered,  i.e.,  (18b),  fourth- 

order  accurate  half-point  differencing  gives 


h. 


h.' 


m 


-h.2  ^ ^ (20a) 

^ [3(bj-bj_^)  '(mj+mj_^)+3(mj-mj_^)  '(b^+b^_^) 

-2  (bm)  j+2  (bm)  ] 


and  (13)  leads  to 


h . 


m . 
3 


-m . , = [b.m.-b.  , , , 

3-1  2 3 3 3-1  3-1  3 3-I- 


, m . , +F  .+F . 


(20b) 


■Y2~  [03m+F)j  -(bm+F)j_j^] 

Once  again  (20a)  equals  (20b)  only  for  constant  b. 

6 . Higher-Order  Extension  of  the  Modified  Finite-Element  Method. 

The  formulation  of  section  3 provides  a consistent  approach  for 
deriving  more  accurate  integral  approximations.  This  can  be  achieved 
by  using  higher-order  polynomials  for  F and  u in  (6) . We  consider 
here  a cubic  spline  polynomial.  Spline  approximations  have  been  used 
earlier  in  differential  collocation  procedures.  The  cubic 

polynomial  approximation  over  the  interval  (*)j_j^,'^^)  can  be  written 
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(1-3) 

as : 

S(T1;3,1)  = SCO)  = hjraj_^t(l-t)-hjrajt^  (1-t) 

(21] 

+Uj_^ (l+2t) (1-t) (3-2t) 

From  (21)  a variety  of  spline  relationships  can  be  derived  ; 

among  these  are 

_ u . -u . , 

M . = ^ (m.  , + 2m.)  - 6 (22j 

3 n.  3-j.  : ^ ^ 

3 

_ U.  --U. 

M . = - ^ (2m. +m.,.)  + 6 1 (221 

J J J"'"-'-  u ^ 

j 

Decrementing  j to  j-1  in  (22)  and  subtracting  the  resulting  equation 
from  (22a),  we  recover  the  two-point  expression  (13), 

h.  h.^ 

Uj-Uj_i  = (mj+mj_^)  - / (23) 

which  forms  the  basis  of  the  fourth-order  extension  of  the  KBS. 

Direct  integration  of  equation  (1)  over  two  adjacent  intervals, 
with  the  polynomial  expression  (21)  for  source  terms,  leads  to 


- li-  + 

From  the  spline  equations  (22a)  and  (22b)  the  following  three-point 
formula  is  obtained: 


XI  . 


11 


r 


Equations  (24)  and  (25)  are  identical  with  the  fourth-order  extension 
of  the  KBS  applied  by  Wornom  (see  equations  (15)  and  (16)).  Wornotn 
eliminates  from  (14)  with  the  non-divergence  form  of  the  equation 
(18b) ; in  the  present  formulation  this  is  unnecessary  as  is  de- 
termined from  the  spline  relationships  (22)  . A small  increase  in  the 
truncation  error  is  introduced  with  (22),  but  the  substitution  of  (18b), 
required  by  Wornom,  is  avoided.  This  can  be  extremely  important  for 
more  complex  equations  and/or  multi-dimensional  flows. 

As  a final  consideration  the  simple  quintic  spline  MFE  is  de- 
veloped. This  has  earlier  been  designated  as  S(T;5,1)  and  leads 

to  a sixth  order  formulation: 


h." 

S(T);5,1)  = S(T)  = Uj_^  + t + t' 


h.' 


+ [lO(Uj-Uj_^)  - hj  (4mj+6mj_^)+  (Mj-3Mj_^)  ] f 


h. 


- [15  (Uj-Uj_^) -hj  (7mj+8mj_^)+  (2Mj-3Mj_j_)  ] t 


+ [6(Uj-Uj_^)-3h^  (mj+mj_^)  + (Mj-Mj_^)]  t^ 

(26) 

With  the  polynomial  (26)  we  integrate  (2a)  over  the  interval 

to  obtain  the  sixth-order  extension  of  the  two-point  Keller 
formula  (13) . We  find 


h.  h/ 


(27) 


h. 


(g;  + g;_i)  + o(h. 


120  ^j-1'  - 

If  we  complete  the  integration  (5)  using  (26)  for  F,  the  sixth-order  MFE 
formula  for  (2a)  becomes 


12 


h . 

h.\  2 ' 2 ' ' , 

- Vi-‘” 

h.3 
+ 120 


In  addition  the  spline  formulas  from  reference  3 are 


7mj^^+8(l+a^)raj+7,.^mj_^=  ^ (Uj^^+(  A ._^)  (29a) 

and 

j 

+ j^(-2m^^^+3(c^-l)mj+2a^m^_^)  (29b) 

Since  F = F(u;Tl),  the  F'  and  F"  terms  in  (28a)  are  at  most  functions 
of  u,m,M.  The  equations  (28),  (29a)  and  (29b)  provide  a (3x3)  block- 

tridiagonal  system  for  (u,m,M)j.  We  recall  that  the  second-order  MFE 
is  a scalar  system  for  u^,  and  the  fourth-order  MFE  is  2x2  for  (u,m)j. 
The  Mj  terms  in  (28,29)  can  be  eliminated  by  using  the  differential 
spline  approximation  to  (1)  in  non-divergence  form.  This  is  the  same 
procedure  used  by  Wornom  to  eliminate  M^  in  the  fourth-order  KBS 
development.  This  was  unnecessary  in  the  fourth-order  MFE  formulation 
as  the  spline  relationships  (22)  were  available  to  evaluate  M^  and 
maintain  a (2x2)  system.  If  a (2x2)  sixth-order  MFE  system  is  pre- 
ferred, it  is  necessary  to  apply  the  non-divergence  M^  elimination* 
however,  as  discussed  earlier  this  is  not  desirable  for  more  complex 
systems  of  equations. 


1 

I 


J 
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7.  Two  dimensional  MFE 


In  this  section,  the  modified  finite  element  method  is  ap- 
plied with  a two-dimensional  equation  of  the  following  type: 

+ (uT)  + (vT)  = (T  +T  ) + q (30) 

t X y Rg  XX  yy  ^ 

If  we  define 


P = q - 


! (31) 


I 

J 


The  above  equation  can  be  written  as: 

= P (32) 

In  this  two-dimensional  case  the  integration  is  over  a rectangu- 
lar grid.  The  method  can  be  used  for  other  element  shapes;  how- 
ever, we  shall  only  consider  a rectangular  element  in  the  present 
investigation.  Bicubic  spline  polynomials  can  be  used  for  curve 
fitting  in  two  dimensions;  however,  in  the  present  analyses, 
partial  splines  are  used  for  discretizing  equation  (32).  This 

amounts  to  a splitting  procedure,  where  the  curve  fit  is  made  in 
the  X direction  and  then  a second  curve  fit  is  made  in  the  y di- 
rection; the  order  is  arbitrary.  This  procedure  has  been  shown 
to  be  equivalent  to  the  use  of  bicubic  splines  (see  reference  11) . 

Integrating  equation  (32)  between  (x^  l'^i+1^'  obtain 

k. 


12 


i-x'  " ^xi-^xi-l’}- 


(33a) 


S . 


1 


-z . 
1 
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where  k.  = x.-x.  , ; c = k..,/k. 

1 1 1-1  X 1+1  1 


S.  = ^ (r-  P.^,+  (l+7  )P.+P.  ,) 
1 2 X 1+1  X 1 1-1 


(33b) 


*TT-  {<Vi+r<Vi 

Integrating  over 


obtain 


2^  "‘"y^^i+l,  j + l'^i-1,  ^^i+1,  j~^i-l,  j) 


+ (F, 


j-l”^i-l,  j-1^  - 


li“  f"y  t*''y>i+l,j+l'*''y’i-l,j  + l"*''y*i+l,j'^<^y*i-l,j  '■ 
k . ^ r 

12  ^^x  i+1,  j+l~  ^'^x^  i,  j+l**^^x^  i+1,  j-l'*"^'^x^  i,  j-1 

{^^x^  i,  j+l~^^x^  i-1,  j+l'^^x^i,  j-l'^^^x^i-l,  j-1 


h. 

= -J.  r 


^ ^-y<2iO+l-2i.j+l^^(l%)  (^iO-^i,j)  ^ (^i,j-l-"iO-in 

-h^/12[,y^{(Sy)i,j,l-(Z^)._^,,-(Sy),^..(Zy).^.} 

*{<Vi,3-'Vi.3-‘Vi.3-l  ^ 'Vio-l}^ 
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7 . 1 Example;  Second-order  MFE  for  Transonic  Small  Disturbance 
Equations 

As  an  example  of  the  two  dimensional  MFE  formulation,  the 

transonic  small  disturbance  equations  are  considered.  In  this 

analysis  only  second-order  accurate  theory  is  discussed;  i.e., 

2 2 

the  and  h^  terms  in  the  pra'i-ious  development  are  neglected. 
The  fourth  order  extension  is  not  considered  here  although  it 
can  be  obtained  from  the  equations  of  section  7.  We  include 


this  example  here  since  the  earlier  formulation 


of  this 


problem  is  only  first-order  accurate  in  the  supersonic  region, 
The  transonic  small  disturbance  equations  are: 


'k  - ( Y+1)  0 T - + * = 0 , 

x-  XX  yy  * 

where  y is  the  ratio  of  specific  heats,  k the  transonic  similar- 
ity parameter  and  ? the  velocity  potential: 


u = ? and  V = 5 
X y 


Rewriting  the  equation  in  conservation  form,  we  obtain 


= 0 , 


This  equation  is  identical  to  equation  (32)  with 


F = ku  - u^ 


G = V 


P = 0 

Since  the  equation  is  of  mixed  type,  the  integration  must  reflect 
the  characteristic  domain  of  dependence  in  the  supersonic  region. 
Therefore  for  subsonic  regions  we  integrate  over  1 ' ^i+1  - ' j-1 j4 l ' 

so  that  (34)  applies  directly.  For  the  siipersonic  domains  we  integrate 
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f 


over  Lx^  2'^i*'  that  (34)  applies  with  i-i-1.  For 

the  full  potential  equations  a similar  procedure  would  apply. 

In  order  to  show  the  second-order  accuracy  for  the  supersonic 
regions,  equation  (33a)  with  i-i-1  becomes,  in  terms  of  (u,v). 


[ku--^2^i  _ ^ 


( ^.x)  . 


i-1 


C " (v  ) . + (l+c  ) (v  ) . 

y 1 X ^ y' i-1 


+ 


^ i-2 


] 


where  - = (ix)./(Ax). 

It  is  easily  shown  that  this  equation  is  second-order  accurate  in 
Ax  even  for  non-uniform  grids.  The  y integration  then  leads  to 
(34).  The  final  expression  for  (u,v)  or  5 involves  the  nine  points 
on  -^j-i'^j+1^'  the  Murman-Cole  first-order 

formula  requires  only  five  points. 

8 . Boundary  Layer  Examples 

The  fourth-order  MFE  method  is  applied  to  the  parabolic  boundary 
layer  equations  written  in  conservation  form.  Comparisons  are  made 
with  second-order  techniques,  the  fourth-order  spline  4 collocation 
method  applied  to  the  non-divergence  differential  form  of  the  equa- 


tions and  the  Wornom  fourth-order  adaption  of  the  KBS.  Specific 
examples  include  (1)  the  similarity  (ordinary  differential) 
equations  governing  the  laminar  flow  over  a flat  plate  and  at  a 
stagnation  point,  (2)  the  quasi-similar  model  equation  for  the 
turbulent  flat  plate  boundary  layer,  and  (3)  the  non-similar 
boundary  layer,  both  laminar  and  turbulent,  described  by  a 
decelerating  linear  external  velocity  field  (Howarth  problem).  In 
almost  all  cases  an  exact  or  very  accurate  numerical  solution  is 
available  for  comparison  purposes. 
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The  governing  boiandary  layer  equations  are 

(d+c)V^)^  + ((f+2?f^)v)^  = V^-0(1-V^)+2?V^  , (35a) 

"n)  = V (.»,?)  = V , (35b) 

where  ? is  the  turbulent  eddy  viscosity  defined  here  by  the 
(9) 

Mrchel  single  layer  model 


.-  = R 


^nl-h 


F = l-exp  — Va"';  a = 26  C(2fR  )^(u^)  =01 

^ e q n ' 


/ = 0.085  T1„(2^R  ) 2rtanh(0.4lT/0.085Tl  )] 

S 0 © ^ 

k ^ 

and  t;  = y(R^/2=-)^  ,.  ? = | u^dx 

'•  o 


2 ? 


du 
€ 

dx 


S(?)  — ~2 — ~ ' ^0  — V = 0.995  and  e denotes  the  boundary 


layer  edge . 

For  the  flat  plate  geometry,  p.=0;  for  the  stagnation  point,  8=1; 
for  the  Howarth  flow,  u = 1-x  and  o = -2 V ( 1-2 for  laminar  flow. 


.-  = 0. 


The  boundary  conditions  are 

f(?.0)  = V(  ?,  0)  =0  and  limV(p,  T])  = 1 . 

T1-.00 


(35c) 
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8.1  Fourth-Order  MFE  Equations 


For  the  boundary  layer  flow  (35)  the  MFE  equations  are  given 
by  (22-25)  with  a » I+0,  b = -(f+2?f,)  and  F = V^-<=(l-v^)  + 2 

for  (35a),  and  a=0,  b=l,  F=V  for  (35b).  The  streamwise 
gradients  are  discretized  by  backward  or  central  differences.  A 
3x3  block-tridiagonal  system  for  (u,M,f)  results. 

The  boundary  conditions  are  given  by  (35c);  in  addition,  the 
boundary  values  for  m or  M are  given  by  (14b),  or  (14b)  and  (22), 
respectively.  An  alternate  boundary  condition,  used  earlier 
for  differential  spline  boundary  layer  calculations  is  given  by 

M(?.0)*=  -o-h2^((26-l)m^(^,0)  + 'm^^(^0))/12  = - 9 h^V^^(  ?,  0) /1 2 . 

where  h2  = TI2  ? denotes  the  grid  point  normal  to  the  surface; 
hj  = is  the  local  mesh  width;  A*'  = 2^  is  the  mesh 

width  along  the  surface. 

8 . 2 Solutions 

(a)  Similar  Laminar  Boundary  Layer 

The  results  for  the  similarity  solutions  ^ (-|^)  =0  in  (35a)_,  for 
flat  plate  (fl=0)  and  stagnation  point  (3=1)  flows  are  given  in 
Tables  1-3  for  a variable  and  uniform  mesh,  with  ten  intervals  be- 


I 


1 

•:'3 


i 


*m(0,?)  can  also  be  obtained  from  a Taylor  series  expansion  at 
T)=0.  For  coarse  grids,  this  can  be  preferable  to  the  spline  re 
lationship  for  m(0,?) 


1 

\ 

i 
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tween  the  surface  and  the  outer  boundary.  For  the  variable  grid 


n = a^c/d+i^o^)  ’ 


(36) 


where  a = 8.26  , b = -0.4  and  a = 24.2538  (1+b  ) 

o o o 


Figure  1 depicts  the  percentage  error  in  wall  shear,  for 
? = 1,  as  a function  of  the  number  of  intervals.  The  velocity 
profiles  for  the  variable  grid  with  o = 0 are  given  in  Table  1 . 
The  KBS4  and  fourth-order  MFE  methods  using  divergence  (conserva- 
tion) form  equations  are  generally  the  most  accurate.  The  KBS4 


method,  which  requires  the  use  of  the  governing  equation  for 


is  slightly  more  accurate  than  MFE;  however,  the  latter  is  less 
time  consuming  and  therefore,  for  equal  accuracy  the  two  procedures 
require  approximately  the  same  CPU  time.  For  more  complex  systems 
of  equations  the  MFE  approach  should  become  more  favorable. 

The  second-order  methods  are  inconclusive.  For  the  flat  plate 
geometry  and  a variable  grid,  the  KBS  is  excellent  and  the  finite- 
difference  poor.  For  f3  = 0,  with  a uniform  grid,  and  for  the  vari- 
able grid  r = 1 solution,  the  opposite  holds;  in  fact,  the  non- 
divergence finite-difference  solutions  are  better  than  the  divergence 
KBS  results.  The  same  is  true  for  the  so-called  Davis  Coupled 


Scheme  (DCS). 


The  solutions  deteriorate,  in  many  cases,  when  non-conservation 
equations  are  considered.  The  KBS2  results  are  particularly  note- 
worthy in  this  respect,  moreover,  the  large  differences  in  the  half- 
point KBS  and  two-point  KBS2  adaption  used  by  Wornotii  are  evident. 
Non-conservation  KBS4  solutions  were  not  available. 

The  non-divergence  spline  4 solutions  are  poor  for  the  variable 
grid  o = 0 case  but  are  significantly  better  for  the  p = 0 uniform 
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grid  and  3=1  variable  grid  examples.  Furthermore,  when  the 
spline  4 method  is  applied  to  the  equations  transformed  by  (36), 
instead  of  a grid  generated  by  (36) , a reasonable  improvement 
results. 

It  should  be  noted  that  the  integral  approach  (PCBS,MFE) 
would  involve  considerably  more  arithmetic  if  (1)  the  transformed 
equations  were  considered,  (2)  the  governing  equations  were  more 
complex  or  extensive  (in  particular,  the  use  of  the  governing 
equation  for  in  KBS4  requires  more  operations) , or  (3)  the 
equations  are  two-dimensional;  as  shown  in  Section  7,  the  complex- 
ity of  the  integrated  equations  is  increased  significantly.  On 
the  other  hand,  the  non-divergence  spline  4 procedure  is  not  modi- 
fied greatly  in  these  cases;  therefore,  the  slight  decrease  in  ac- 
curacy may  be  offset  by  its  simplicity.  Finally,  the  spline  4 
procedure  can  be  adapted  in  order  to  obtain  conservation  solutions 
in  differential  form;  in  references  (2,3)  significant  improvements 
in  accuracy  were  noted.  This  requires  several  additional  curve 
fits,  but  does  not  appear  to  be  as  complicated  as  the  two-dimensional 
KBS  or  MFE  procedures  for  simple  rectangular  elements.  For  other 
element  shapes,  the  integral  finite-element  approach  might  be  more 
preferable . 

In  summary,  it  can  be  concluded  that  the  accuracy  of  the  results, 
of  any  of  the  fourth-order  methods  described  here, is  very  dependent 
on  the  combination  of  pressure  gradient  parameter  0 and  the  choice 
of  grid. 
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(b)  Quasi-Similar  Turbulent  Boundary  Layer 

For  these  calculations  we  assume  that  ^ = 0 in  (35a) , even 
though  the  solutions  do  depend  on  (Rg-).  This  is  a fair  approxi- 
mation and  allows  us  to  consider  an  ordinary  differential  equation 
for  the  turbulent  boundary  layer.  In  subsequent  examples,  the 
full  non-similar  equations  will  be  considered. 

The  mesh  is  given  by  the  grid  generated  from  (36) , with 

1 = -109,  b = 0.05  and  a = 60(l+b  ) '^\  The  percentage  error  in 
o o o 

shear  stress  is  depicted  on  Figures  2 and  3,  for  S = 0 and  p = 0.5, 
respectively.  For  these  cases  in  the  turbulent  eddy  viscosity 
model  is  given  as  24.5  and  15.8,^^^  respectively.  The  results  for  the 
KBS4  and  MFE  are  good,  although  less  than  fourth-order  accurate. 

This  appears  to  be  caused  by  the  error  in  the  numerical  evaluation 
of  the  eddy  viscosity  =:.  For  fixed  e,  fourth-order  accuracy  should 
be  recovered. 

The  non-divergence  form  solutions,  obtained  with  the  grid 
transformation  are  reasonable,  but  somewhat  erratic.  The  second- 
order  finite-difference  solutions  are,  surprisingly,  more  accurate 
than  the  divergence  KBS2  results.  The  non-divergence  KBS2  solutions 
are  rather  poor. 

(c)  Non-Similar  Solutions 

For  the  non-similar  boundary  layer  calculations  a,b  and  F in 


(35)  are  defined  in  section  8.1.  Three  cases  have  been  considered 


here:  (1)  turbulent  flow  with  B = 0.5;  (2)  laminar  Howarth  problem, 

and  (3)  turbulent  Howarth  problem.  In  the  first  case  the  quasi- 
similar solution  has  previously  been  discussed.  The  effects  of  the 
non-similar  terms  will  be  examined.  For  case  (2)  the  separation 


1 point  has  been  critically  examined  in  many  investigations;  it  has 

I 

1 


T 
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been  found  that  separation  occurs  at  x = 0.1198.  In  this  section, 

separation  point  calculations  for  course  grids  are  compared  with 

(12) 

earlier  solutions  obtained  with  the  non-divergence  spline  4 

system  of  equations,  as  well  as  second-order  methods.  Finally, 
the  suppression  of  separation  for  the  turbulent  boundary  layer  in 
a decreasing  linear  external  flow  is  described.  To  the  authors' 
knowledge  solutions  for  this  flow  have  not  been  published.  This 
case  is  considered  here,  simply  to  demonstrate  the  applicability 
of  the  MFE  method  to  turbulent  boundary  layers  in  adverse  pressure 
gradients . 

The  solutions  are  shown  in  Tables  4-6.  For  constant  p = 0.5 
the  grid  defined  by  (36)  with  10  intervals  in  11  is  specified; 

■'^^=15. 8 in  the  eddy  viscosity  model.  This  case  has  previously 
been  treated  by  Wornom  using  the  KBS4  procedure.  The  effect 

of  the  marching  increment  -i?  and  comparisons  with  the  quasi-similar 
solution  at  ^=1.88x10^  are  given  in  Table  4.  For  Wornom 's  grid 
the  are  given  as:  0,  86xl0~^,  43xl0~^,  86x10”'^,  43xl0~^,  0.17, 

0.26,  0.34,  0.51,  0.68,  0.86,  1.03,  1.28,  1.50,  1.71  and  1.88  for 
i=l,...,17.  There  is  a small  effect  of  the  ? spacing,  but  more 
important  is  the  assumption  that  T|^=15.8  throughout.  This  should 
only  be  true  at  R^'=1.89{10  , and  therefore  the  final  solutions  have 
an  inherent  error  associated  with  the  inaccurate  71^  value  used  for 
?<1.88.  A more  accurate  solution*  would  require  the  precise  evalua- 
tion of  at  each  ? location;  however,  with  a crude  grid  this  is 

*In  a more  exact  calculation  with  60  intervals  it  was  found  that  at 

R =430,000,  T)  «10. 
e e 
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not  possible,  and  this  leads  us  to  believe  that  for  the  turbulent 
boundary  layer,  the  error  estimates  obtained  with  the  model  quasi- 
similar problem  are  somewhat  artificial. 

For  the  laminar  Howarth  problem  a very  coarse  uniform  grid, 

with  h=1.0  provides  reasonably  accurate  estimates  of  the  separation 

point,  see  Table  5.  The  MFE  results  are  somewhat  less  accurate  than 

{ 12) 

those  obtaxned  with  the  non-divergence  spline  4 method  but 

both  are  significantly  better  than  the  second-order  solutions. 

Some  calculations  were  also  made  with  the  variable  grid  considered 
for  the  similarity  solutions.  As  was  found  with  the  similarity 
problems,  the  non-divergence  spline  4 results  become  less  accurate 
than  those  obtained  with  the  uniform  grid.  Surprisingly,  the  MFE 
solutions  also  deteriorate;  and,  for  a reason  that  is  not  clear, 
these  calculations  require  orders  of  magnitude  more  iterations  for 
convergence  than  do  the  spline  4 calculations  or  the  MFE  solutions 
with  the  uniform  grid.  It  would  appear  that  the  MFE  calculations 
are  somewhat  more  sensitive  to  the  grid  generated  by  (36)  when  non- 
similar flows  with  adverse  pressure  gradients  are  considered.  The 
similarity  results  using  MFE  and  (36)  were  extremely  accurate  and 
rapidly  convergent. 

For  the  turbulent  Howarth  problem,  separation  does  not  occur 
for  »<0.31  or  Rg?<580,000;  the  value  in  the  eddy  viscosity  model 
is  evaluated  at  eachf  location  and  the  transformation  (36)  with  60 
intervals  is  specified  in  order  to  accurately  estimate  The  MFE 

results,  on  Table  6,  appear  to  be  reasonable,  although  the  accuracy 
is  difficult  to  assess.  The  “ grid  is  rather  crude,  a«:0.005,  and 
comparisons  with  other  procedures  were  not  possible. 
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In  summary,  the  accuracy  of  the  various  methods  described 
herein  appears  to  be  dependent  on  a combination  of  (1)  the  prob- 
lem, i.e.  laminar  or  turbulent,  constant  g or  variable  g,  (2), 
the  choice  of  the  grid  and  (3)  non-divergence  or  divergence  form 
of  the  equations.  Generally,  the  fourth-order  methods  lead  to 
improvements  over  second-order  techniques  but  the  relative  accuracy 
of  the  fourth-order  procedure  is  dependent  on  the  factors  just 
mentioned . 


25 


References 


1.  Rubin,  S.G.  and  Khosla,  P.K.:  "Higher  order  numerical  solu- 

tions using  cubic  splines".  AIAA  Journal  Vol.  14  No.  7 1976 
pp.  851-858,  also  NASA  CR2653. 

2.  Rubin,  S.G.  and  Khosla,  P.K.:  "Numerical  methods  based  on 

polynomial  spline  interpolation".  Proceedings  of  Fifth  Int'l 
conference  on  Numerical  Methods  in  Fluid  Dynamics,  Enschede, 

The  Netherlands,  Springer-Verlag,  Heidelberg,  1976  pp.  370-377. 

3.  Rubin,  S.G.  and  Khosla,  P.K.:  "Polynomial  interpolation  methods 

for  viscous  flow  calculations".  To  appear  in  the  Journal  of 
Computational  Physics.  Also  NASA  CR2735. 

4.  Keller,  H.B.:  "A  new  difference  scheme  for  parabolic  problems". 

Numerical  solution  of  partial  differential  equations",  Vol.  II 
edited  by  J.  Bramble,  Academic  Press,  New  York  1970. 

5.  Wornom,  Stephen  F.:  "A  critical  study  of  higher-order  numerical 

methods  for  solving  the  boundary-layer  equations".  AIAA  3rd 
Computational  Fluid  Dynamics  Conference,  Albuquerque,  New  Mexico, 
June  1977.  Also  NASA  TMX-74000. 

6.  Ackerberg,  R.C.  and  Phillips,  J.H.;  "The  unsteady  laminar 
boundary  layer  on  a semi-infinite  flat  plate  due  to  small  fluctu- 
ations in  the  magnitude  of  the  frees tream  velocity".  J.  Fluid 
Mech.  Vol.  51,  pp.  137-157  (1972). 

7.  Blottner,  F.G.;  "Investigation  of  some  finite-difference  tech- 
niques for  solving  the  boundary  layer  equations".  Computer  methods 
in  Applied  Mechanics  and  Engineering,  Vol.  6,  pp.  1-30,  1975. 


I 


I 


26 


8.  Keller,  H.B.;  "Numerical  solution  of  boundary  value  problems 

for  ordinary  differential  equations":  survey  and  some  recent 

results  on  difference  methods.  Numerical  solutions  of  boundary 
value  problems  for  ordinary  differential  equations,  A.  K.  Aziz, 
ed..  Academic  Press,  New  York,  pp.  27-88  1975. 

9.  Michel,  R,;  Qudmard,  C.;  and  Durant,  R.;  "Hypothesis  on  the 
Mixing  length  and  application  to  the  calculation  of  the  turbulent 
boiindary  layers".  Proceedings  of  Computation  of  Turbulent 
Boundary  layers  - 1968  AFOSR-IFP-Stanford  Conference,  Vol.  1, 

pp.  195-212. 

10.  Rubin,  S.G.,  Khosla,  P.K.  and  Rivera,  S.:  "Turbulent  boundary- 

layer  studies  using  polynomial  interpolation".  Presented  at 
Symposium  on  Turbulent  Shear  Flows  at  Penn  State  University, 

Penna . , April  18-20,  1977.  (To  be  published). 

11.  Ahlberg,  J.,  Nilson,  E.  and  Walsh,  J.:  "Theory  of  splines  and 

their  applications".  Academic  Press,  New  York,  1967. 

12.  Riibin,  S.G.  and  Rivera,  S.:  "Separating  and  turbulent  boundary 

layer  calculations  using  polynomial  interpolation."  Polytechnic 
Institute  of  New  York  Report  No.  POLY  M/AE  77-4.  February  1977. 

13.  Murman,  Earl  M.  and  Cole,  Julian  D.:  "Calculation  of  Plane  Steady 

Transonic  Flows".  AIAA  J.  9,  1,  (1971)  pp.  114-121. 


27 


Table  1:  Laminar  Flat  Plate  ( :=0) : Variable  Grid  to  24.2538 

(10  intervals) 


(a)  Surface  Shear  Stress 

Conservation  Form  Solutions 


Exact 

MFE 

KBS4 

KBS2 

F.D. 

0.469600 

0.470362 

0.469835 

0.468759 

0.380908 

Non-Conservation  Form  Solutions 


SPLINE 4 

TRANSFORMED 

PLANE 

KBS 

KBS2 

DCS 

SPLINE4 

0.438678 

0.476971 

0,494798 

0.421500 

0.488898 

(b) 

Velocity  Profiles 

Conservation  Form 

Solutions 

Grid 

Exact 

KBS4 

MFE 

KBS  2 

0.0000 

.00000. 

0,00000 

0.0000000 

0 

.00000 

0.0500 

.02347 

.02348 

,0235044 

.02343 

0.1421 

.06671 

.06674 

.0667994 

.06659 

0.3076 

.14437 

.14444 

,1446144 

.14412 

0.6023 

.28164 

.28178 

,2819718 

,28113 

1.1266 

.51474 

.51498 

.5157463 

,51364 

2.0651 

.83286 

.83293 

.8350900 

.83263 

3.7657 

.99552 

.99616 

.9935382 

1 

.01397 

6.9004 

1.00000 

.99947 

1.0010271 

.99643 

12.8087 

1.00000 

.99983 

.9995595 

1 

,00320 

24.2538 

1.00000 

1.00000 

1.0000000 

1 

.00000 

Non-Conservation  Form  Solutions 

SPLINE4 

TRANSFORMED 

KBS 

KBS2 

DCS 

PLANE 

SPLINE4 

.00000000 

0.0000 

0.0000 

0.00000 

0 

.0000 

.02192272 

.0246 

.0247 

.02106 

.0244 

.06231277 

,0699 

.0703 

,05987 

.0694 

.13487369 

.1513 

.1521 

.12956 

.1503 

.26319919 

.2949 

.2963 

.25264 

.2936 

.48176760 

.5396 

.5382 

.46098 

.5393 

.78493782 

.8605 

.8544 

.75179 

.8929 

.93994754 

1.0187 

1.0156 

.97564 

1 

.0820 

.96388488 

,9999 

.9960 

1.00253 

.9969 

.95144079 

1.0259 

1.0036 

.99166 

1 

.0687 

.00000000 

1.0000 

1.0000 

1.00000 

1 

.0000 

Table  2: 

Laminar  Flat  Plate 

(a=0);  Uniform  Grid  (h=l) 

Exact 

(NC)  S4 

(NC)  KBS2 

(NC) FD 

(0.24%) 

(6.15%) 

(1.61%) 

0.469600 

0.470730 

0.440743 

0.47718 

jMFE2( 

MFE4 

KBS4 

|kbs2| 

(6.9%) 

(0.003%) 

Unavilable 

0.437207 

0.469614 

Table  3: 

Laminar  Stagnation  Point  (?=1):  Variable  Grid 

24.2538  (10  intervals) 

Exact 

(NC)F.D. 

(NC)S4 

(NC)DCS 

( 1 . 17%) 

(0.141%) 

(1.53%) 

1.23259 

1.24826 

1.23084 

1.25140 

F.D. 

DCS 

MFE 

KBS4 

(0.226%) 

(2.39%) 

(0.042%) 

(0.016%) 

1.23538 

(1.26205) 

1.23311 

(1.2328) 
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Table  4 


Turbulent  Non-Similar  Solutions  (o=0.5,  R =1.88x10  , = 

e z e 

(a)  X 10^ 


MODEL 

MFE  (A— 0.005) 

MFE  (WOPUOM  GRID 

3.5473 

3.9182 

3.9344 

(b)  VELOCITY 

PROFILE 

GRID 

MODEL 

MFE  ( A'='=0.005) 

MFE(Vv’ORNOM  GRID) 

0.0000 

0.0000 

0.0000 

0.0000 

0.0500 

0.1715 

0.1828 

0.1836 

0.1421 

0.3742 

0.3935 

0.3948 

0.3076 

0.5162 

0.5387 

0.5401 

0.6023 

0.6052 

0.6296 

0.6312 

1.1266 

0.6716 

0.6980 

0.6997 

2.0651 

0.7342 

0.7624 

0.7643 

3.7657 

0.7971 

0.8269 

0.8288 

6.9004 

0.8801 

0.9110 

0.9130 

12.8087 

0.9636 

0.9810 

0.9814 

24.2538 

1.0000 

1.0000 

1.0000 

Table  5: 

Laminar  Howarth  Problem:  x 

separation 

(NC) 

MFE  2 

Exact 

Spline4 

MFE  4 

KBS2  F.D. 

( 

} 

Uniform  Grid 
h=1.0 

0.1196 

0.1225 

0.1374  0.1495 

Variable  Grid 

0.1096 

0.1328 

0.1515  

0.1198 

7 =24.2538 

max 

30 

L 


15.8) 


Table  6:  Turbulent  Howarth  Problem  (Re=l ,88x10^) 

(a)  Skin  Friction  Coefficent 


0.0410 

0.0710 

0.1010 

0.1310 

0.1610 

CjXlO^ 

6.1360 

5.2861 

4.7318 

L- — 

4.2847 

3.8721 

=• 

0.1910 

0.2210 

0.2510 

0.2810 

0.3110 

CjXlO^ 

3.4849 

3.0593 

2.5985 

2.0266 

1.2409 

(b)  Velocity  Profile  R^_=584680 

Tl^=23 .3468 


Grid  Location 

0.0 

0.'?3f'7B?71  0-02 
0.117SS137O-01 
0.  iqi05691O-.j  I 
0.2^1B10H90-01 
0. 38'^62^23l)-ul 
0,50^5<i-658n-ol 
0.64^fl4642D-i)l 
0.^09089420-01 
0.9U813913Q-01 
0.12112005D  00 
0.14S73572D  00 
O.17401132Q  00 
0.2n^2‘t.tt000  00 
O.24?90282D  00 
0.2S473<*430  00 
0.332269390  00 
0.3841291'iO  00 
0.447133910  00 
0.516140790  00 
0.59410363D  00 
0.6«20H3170  00 
0.7P125656D  00 
0.«9?94008n  00 
0.101858300  01 
0.115980330  01 
0.131339410  01 
0.14Q634510  01 
0.16q586260  01 
0.191939290  01 
0.216964680  01 


Velocity  Profile 
O.T) 

0 , 3f'25ri40  ID  - Or:' 

0 .‘^0  00332‘*0 - 02 
0.1  329893  ^O-O 1 
0 . 19-,59r,l  ID-ul 
0,27U9309qD-0l 

0,360861rt4Q-(i  1 

0 .464086220-0 1 
0.595739r,30-01 
0.7<*7294S7r)-Ol 
0.926221 700-01 
0.11J51693D  00 
0.137461670  00 
0.164136980  00 
0.192821600  00 
0. 222545020  00 
0.252352710  00 
0.281489610  00 
0,309470670  00 
0.336037310  00 
0.361112640  on 
0 . 384720060  00 
0.406965840  00 
0.427972810  00 
0.4+7905800  00 
0.4669o31«0  00 
0.485142360  00 
0.502746040  00 
0.519895290  00 
0.536646040  GO 
0.553309500  00 


Grid  Location 
0.244962780  01 
0.276266211  01 
0.3112432^0  01 
0.350301470  01 

0.393891910  111 

0 . 4 4 8 5 1 3 4 2 0 0 1 

0.496717680  01 

0.557114550  Ul 
0.624378050  01 
0.699252940  01 

0.7«256183D  01 
0.87P21320D  01 
0.9782inoOn  01 
0.109265930  02 
0.121978260  02 
0.136092770  02 
0.151758110  02 
0.169138180  02 
0.18»41370f)  02 
0.209753840  02 
0.233468090  02 

0.259708280  02 
0.2««77075O  02 
0.320948860  02 
0.3'^656556n  02 
0.395976410  02 
0,430572700  02 
0.487785020  02 
0.541087100  02 
0.600000000  0? 


Velocity  Profile 

O.564‘^2'''64  0 uO 

0 , 5 8 6 4 5 9 4 — u I ' 

0 . 6 0 32  0 30  (I  n 

0 . 6 0 3 ’•<  - 1 3 L)  " 

0 . 6 3 7 9 I 3 - u I ■ ( 1 

0 . 6'^ 49 2 2450  f,  0 

0.64475977D  uO 
0 . 7 1 'D  38652D  no 
0,737-265t)O  no 
0,760422580  on 
0.788207050  on 
0.81Q93  + U6L'  00 
0.838680900  00 

0.868831940  ,)0 

0.8466'’0580  go 
0.9?8;?1O3mD  go 
0.9535-4110  00 
0 . 4 7 8 6 1 H 1 4 0 0 0 
0,99+515-20  on 

0 , 9 44  9 3[;  0 0 

0.1  (JO  00  2 320  0 1 
0.9499476  + 1.1  uO 
0 . 1 (J  0 0 0 0 0 2 JJ  0 1 
0. 4449994 -D  00 
0 . 1 0 u 0 (^  0 0 0 0 0 1 
0.10(30  0 0000  0 1 
O.IOOOOOOUI;  0 1 
0 . 10  0 0 0 0 u u in  n 1 
O.IOOGOOOOO  01 
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PERCENT  ERROR  IN  WALL  SHEAR 
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FIG.1  PERCENT  ERROR  IN  THE  WALL  SHEAR: 
STAGNATION  POINT  FLOW  (i9-l) 
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PERCENT  ERROR  IN  WALL  SHEAR 
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FIG.  2 PERCENT  ERROR  IN  THE  WALL  SHEAR 
FOR  MODEL  PROBLEM 


PERCENT  ERROR  IN  WALL  SHEAR 
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FIG.3  PERCENT  ERROR  IN  THE  WALL  SHEAR 
IN  TURBULENT  MODEL  EQUATION  (/3«0.5) 
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